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INTEGRATED ANALYSIS AND DESIGN OF THICK COMPOSITE STRUCTURES 
FOR OPTIMAL PASSIVE DAMPING CHARACTERISTICS 


D.A. Saravanos 
Ohio Aerospace Institute 
Brook Park, Ohio 44142 


ABSTRACT 


The development of novel composite mechanics for the analysis of damping in composite 
laminates and structures, and the more significant results of this effort are summarized. 

Laminate mechanics based on piecewise continuous in-plane displacement fields are 
described that can represent both intralaminar stresses and interlaminar shear stresses 
and the associated effects on the stiffness and damping characteristics of a composite 
laminate. Among other features, the mechanics can accurately model the static and 
damped dynamic response of either thin or thick composite laminates, as well as, 
specialty laminates with embedded compliant damping layers. 

The discrete laminate damping theory is further incorporated into structural analysis 
methods. In this context, an exact semi-analytical method for the simulation of the 
damped dynamic response of composite plates was developed. A finite element based 
method and a specialty four-node plate element was also developed for the analysis of 
composite structures of variable shape and boundary conditions. Numerous evaluations 
and applications demonstrate the quality and superiority of the mechanics in predicting 
the damped dynamic characteristics of composite structures. 

Finally, additional development was focused towards the development of optimal tailoring 
methods for the design of thick composite structures based on the developed analytical 
capability. Applications on composite plates illustrated the influence of composite 
mechanics in the optimal design of composites and the potential for significant deviations 
in the resultant designs when more simplified (classical) laminate theories are used. 
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1. SUMMARY OF ORIGINALLY PROPOSED RESEARCH 


The originally proposed research plan addressed the need to develop novel damping 
mechanics for the accurate analysis and design of thick composite structures. The 
proposed plan emphasized the need of damping mechanics applicable not only on thin 
composite structures, but also on composite structures of significant thickness; structures 
with embedded compliant interlaminar layers; and structures at elevated temperatures. 

In order to realize the full benefits of damping in the analysis and design of advanced 
composite structures, the following technical objectives were defined: 

(1) Development of admissible integrated damping mechanics for the accurate 
simulation of damping in thick composite laminates and/or laminates with interla- 
minar damping layers. 

(2) Development of novel mixed structural mechanics enabling the accurate 
computational simulation of damping in thick composite structures, at room or 
elevated temperature conditions. 

(3) Development of novel tailoring methods for the design of thick composite 
structures. The methods were proposed to integrate the damping mechanics with 
suitable optimization methodologies, and address current structural requirements 
via multiple design criteria. 


2. SUMMARY OF SIGNIFICANT ACCOMPLISHMENTS 
2.1. Research and Technology Developments 

The three proposed milestones of the proposed research plan were accomplished as 
follows: 

(1) Thick-Laminate Damping Mechanics: The critical issue of accurate representation of 
displacement and stress variations through the composite plies was addressed by 
developing rigorous discrete laminate damping mechanics incorporating piecewise contin- 
uous displacement fields (refer to Section 3; Fig. 3.1). The effects of temperature and 
moisture were incorporated. The mechanics were evaluated on thick laminates, laminates 
with interleave layers, and laminates at higher temperatures (for more details refer to 
Sections 4 and 5; Figs. 4.1, 4.2, 5.3-S.6). In addition to superior damping predictions, 
more accurate evaluations of laminate displacements and strains were obtained, which 
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significantly enhanced the in-house modelling capacity for thick laminates (see also, code 
ICAN/DLDT, Section 2.3). 

(2) Computational Simulation of Structural Damping: The developed thick-laminate 
damping theory was integrated into a finite element discretization scheme. A specialty 
plate finite element was developed, enabling the synthesis of damping, stiffness, and mass 
matrices (refer to Section 3). 

The element formulation allows application in two different wa}^ specifically: (1) as a 
generalized plate element with 4 nodes at the middle surface of the composite; and (2) 
as a generalized "thick-shell” type element with nodes at the upper, lower, and optionally 
middle surfaces of the laminate. Evaluations of this novel finite element, has shown great 
advantages in the simulation of coupled thermomechanical response, not only in thick 
composite structures (refer to Section 4; Figs. 4.1, 4.4-4.8), but also, in composites with 
compliant interply layers (refer tO' Section 6; Figs 6. 1-6.8) and thin composite structures 
with significant interlaminar inhomogeneity (see Fig. 4.2). The element is presently 
integrated in the STAT (Structural Tailoring of Advanced Turboprops) in-house code. 

(3) Optimal Design of Thick Composite Structures: A design method for the tailoring of 
composite components of general thickness, lamination, geometry and loading/boundary 
conditions was developed. The method is integrated in the STAT (Structural Tailoring of 
Advanced Turboprops) in-house code (code STAT/DLDT, Sec. 2.3). In this manner, the 
tailoring capability of the code is extended to cover the area of thick composite 
structures. The design methodology is consistent with previous developments addressing 
current challenging structural requirements expressed as multiple design criteria. A multi- 
objective formulation is available and candidate design variables may include parameters 
from the multiple structural scales of a composite structure. 

An additional smaller scale method was also developed for tafloring thick composite 
plates. The method is encoded into a stand-alone code (code THICOPT, Section 2.3) 
which can provide quick, inexpensive optimal designs for simply-supported composite 
plates. Applications and comparisons with similar design methods based on thin 
composite mechanics demonstrated significant advantages of the new method (refer to 
Section 7; also see Fig. 7.1). 

In both cases mentioned above, the optimization is numerically performed with direct 
non-linear programming algorithms, that is, algorithms performing a direct search in the 
design domain based on first order sensitivity. 
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2.2. Summary of Benefits and Results 
The following significant results were attained: 


(1) Development of functional in-house software. Specifically the following codes 
were developed: 

ICAN/DLDT: An integrated code for the dynamic analysis of thick 

composite laminates and plates incorporating the developed 
discrete laminate damping theory (DLDT). An addendum to 
ICAN’s Usser’r Manual is under preparation 


STAT/DLDT: A structural analysis/tailoring code for composite propulsion 

components. The code now includes the developed finite 
element as previously described 

THICOPT: (Thick Composite Plate Tailoring). A versatile, 

computationally inexpensive code for optimal design 
applications in thick composite plates. 


(2) Evaluation and demonstration of the quality of the developed novel damping 
mechanics on select thick composite structures. The significant advantages of the 
developed laminate and structural damping mechanics were quantified for 
composites structures of: (i) higher thickness; (ii) alternating angle-ply laminates; 
and (iii) laminates with compliant interlaminar layers. 

(3) Initial evaluations of the developed optimal design methodology on typical 
composite plates were performed. The results demonstrated significant benefits in 
the design of composite structures. 

(4) Numerous technical publications (see subsection 2.4) describing the developed 
methodologies, and enumerating evaluations, experimental correlations and 
conclusions of technical significance. 


23. Resultant Publications 

Saravanos, D. A. and Pereira J. M., 'Effects of Interply Damping Layers on the Dynamic 
response of Composite Plates," AIAA Journal, in press, Nov. 1992. 
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Saravanos, D. A., “Simulation of Passive Damping in Thick Composite Structures,” 
AIAA J.. to appear. 

Saravanos, D. A., “Damped Dynamic Characteristics of Thick Composite Laminates and 
Plates,” ASME Journal of Applied Mechanics, to appear. 

Saravanos, D. A. and Chamis C. C., "Damped Dynamic Characteristics of Specialty 
Composite Materials and Structures," Proceedings. Symposium on Vibroacoustic 
CharacteriTatinn of Materials and Structures. ASME Winter Annual Meeting, Anaheim, 
California, Nov. 8-13, 1992 

Saravanos, D, A. "Analysis and Optimal Design of Thick Composite Structures with 
Passive Damping Considerations," Proceedings. 4th AIAA/AI R FORCE/NASA/OAI 
Symposium on Multidisciplinary Analysis and Optimization. Sept. 21-23, 1992, Qeveland, 
Ohio. 
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3. THEORETICAL BACKGROUND 


3.1. Introduction 

The significance of passive damping in improving the dynamic performance of flexible 
structures requiring tight vibration control, high fatigue endurance, and accurate position- 
ing of devices and sensors has been recently emphasized. In addition, low levels of 
damping are also considered beneficial in actively controlled structures, as they typically 
improve the robustness of the controller. As a result, the prediction and tailoring of the 
passive damping capacity of polymer matrix composite structures is receiving current 
attention. 

Various damping mechanics theories for unidirectional composites, laminates and 
structures have been reported and a representative selection is mentioned here^ “. 
However, most of the previously mentioned work is limited to the classical laminate 
theory assumptions which neglect interlaminar shear effects, hence, this work is suitable 
only for thin laminates. Moreover, interlaminar shear effects are more critical in 
laminated composite structures because of their through-the-thickness inhomogeneity and 
anisotropy, than in homogeneous isotropic materials. Therefore, interlaminar shear 
damping is expected to be significant in composites, particularly in laminates of thicker 
sections. Yet, limited research has been reported in this direction^^. 

This paper presents recent developments on discrete damping mechanics enabling more 
accurate prediction of damping properties in thick composite structures. The modelling 
of damping at the laminate level is based on a novel discrete-layer laminate damping 
theory (DLDT) which assumes a discrete, yet, piecewise continuous displacement field 
through the laminate. The DLDT effectively captures both interply and intraply stresses. 
The overall damping capacity of the laminate includes contributions from extension, 
flexure-shear, and various material coupling from possible asyimnetries that may exist due 
to the heterogeneous nature of the composite laminate. 

Semi-analytical methods were first developed for the prediction of modal damping and 
other dynamic characteristics of thick specialty composite plates^^ This method was exact 
and computationally inexpensive, but its applications were limited. In order to simulate 
damping in composite structures of general thickness, shape, lamination, and boundary 
conditions, generalized structural mechanics based on finite element discretization were 
developed. This approach is described in the following sections. A specialty 4-node 
element incorporating the aforementioned discrete-layer laminate damping theory is 
described. Generalized damping, stiffness and mass element matrices are formulated. 
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Therefore, in addition to improved modal damping predictions, more accurate 
calculations of the dynamic response and stresses may be also obtained. 

Evaluations of the damping mechanics are presented. The validations include either 
direct comparisons with known exact solutions for composite plates, or comparisons 
between the DLDT and classical damping theory. The significance of interlaminar shear 
damping contributions is also quantified. Finally, the effects of thickness and lamination 
on the damped dynamic characteristics of composite plates and beams are investigated. 


3 . 2 . Composite Mechanics 

This section briefly reviews the synthesis of damping for on-axis composites (damping 
along the material axes), off-axis composites (unidirectional composites loaded at an 
angle), and composite laminates of general stacking sequence. The measures of damping 
used in this paper will be primarily the loss factor rj> ^ind secondarily the specific damping 
capacity (SDC) »|r (r? = i|r/ 27 r). 

Unidirectional Composites 

For a unidirectional composite loaded along the material axes (see Fig. 3 . 1 a), closed-form 
expressions have been developed for the synthesis of elastic and dissipative properties^. 

In summary, five independent elastic parameters completely characterize the stiffness of 
a unidirectional composite (orthotropic but transversely isotropic material). An additional 
four independent damping loss factors characterize the composite damping, that is, 
longitudinal loss factor rjni (direction 11 ), transverse in-plane damping r)^ (direction 22 ), 
transverse through-the-thickness damping 77133=^122 (direction 33 ), in-plane shear damping 
77,66 (direction 12 ), interlaminar shear damping 77,44 (direction 23 ), and interlaminar shear 
damping 77,55= (direction 13 ). Both damping and elastic properties are explicitly 
related to the fiber and matrix properties, and to the fiber volume ratio (FVR). 

For the case of off-axis composites, ie. composites loaded at an angle 0 , a 6 by 6 
damping matrix [77J best describes the damping of the composite. More details about this 
off-axis damping matrix are provided in the Appendix II and refs. 7 , 13 . Off-axis loading 
affects the overall damping capacity of the composite in two distinct ways that are 
uncommon to isotropic materials: by altering the dissipative capability of the ply directly 
associated to normal and shear strains; and by inducing damping coupling between 
normal and shear strains. 
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Thick Laminates 


To model the damping of thick composite laminates, a discrete layer laminate damping 
theory (DLDT) incorporating a piecewise continuous displacement field through-the- 
thickness is developed. The kinematic assumptions for the laminate theory are 
schematically shown in Fig 3.1b. Discrete layer or generalized laminate theories have 
been proposed (see refs. 14, 15) for the more accurate calculation of stresses in thick 
laminates. However, with additional developments presented herein, the DLDT combines 
the potential for accurate damping predictions in composite laminates while maintaining 
generality and elegance. The assumed displacement field has the form, 

u 

v(x,y,z/)= v®(W) +v(w^) (1) 


where superscript o represents the uniform through-the-thickness midplane displacement, 
and the bar indicates the through-the-thickness variation in the displacement field. 
Assuming that the displacements are separable functions of z, the previous equations 
take the form [15]: 


u{xy^4)= 






/-I 




v°(W)+ 
w °(x^^) 


N 




( 2 ) 


where, u* and v* are displacements, along the x and y directions respectively, at the 
interfaces between composite plies or sublaminates (group of plies). P(z) are 
interpolation functions. In this manner, the assumed in-plane displacement field is 
general, in that it may represent extensional, flexural, shear, and coupled deformations as 
well as interlaminar shear strains through the thickness of the laminate. 


The engineering strains {e^} in each composite ply are directly derived fi"om eqs. (2) 
The normal strain is zero as a result of the assumed constant through-the-thickness 
deflection w. The midplane strains are given by, 
and the generalized strains are. 

The comma in the subscripts indicates differentiation. 
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;-i 

e^= 0 


O . ^ 

«c7=«^ 


o c> 


o o 


o o 


o o o 
€ ,=W +V 


<r «; 

4-''' 


4' K 

4» «'■ 


4-«; 


The dissipated strain energy per unit area of the laminate AS^ is, 


AS^ = l/lQ 2n e/[QMe^ 


( 3 ) 


(4) 


(5) 


( 6 ) 


where [QJ is the off-axis composite stiffness matrix. Combination of eqs. (6, 2-5) 
ultimately provides the dissipated strain energy per unit area. 


N 


N N 


AS, = hn(ef\AX * * E E 

2 7«i y»i w»i 


(7) 


where [Aj] is the extensional laminate damping matrix, which includes out-of-plane shear 
terms. The generalized coupling damping matrices [B^j] and flexural/shear matrices [D</“] 
are new. The matrix expressions in eq. (7) are given by 

Similarly, the maximum laminate strain energy per unit area Sl is by definition. 


Sr= 


ha 


J -*/2 




(9) 


Combination of eqs. (3-5, 9) provides the maximum strain energy as a separable form of 
material properties and strains, 

where, [A] is the extensional laminate stiffness matrix with additional out-of-plane shear 
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i,n =1,2,6 

E [ ^ ([ecL[n.],)^i(z>fe 

//I =4,5 


i,n =1,2,6 

»-i JVi 

z,n =4,5 

N N N 


S^= l/2(6fp4K - 2cfE [B'V. » E E 


y«l y«l m«l 


( 8 ) 


( 10 ) 


terms, {&] are the generalized coupling stiffness matrices, and the flexural/shear 
matrices^^. These matrix expressions are given by: 


t4]= E [QcV^ 

k~l >k-l 


*.l J*k-1 

i/i =1,2,6 

ifc-1 J^k-l 

i/i =4,5 

^1 

E f** ([2cL)»/'W"(^)* 

H-l J**-! 

i/i =1,2,6 

k-i 

i/i=4,5 


Considering eq. (2), the laminate kinetic energy per unit area takes the form: 


( 12 ) 
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where the generalized laminate mass and inertia matrices are defined as, 




S diag{p,)dz [Bil 

fc-i 

I. 

diag{p,)F\z)F”^{z)dz 

k-l 


^1 


E 


[*'• diag((,,)F\z)dz 
J*k-1 


( 13 ) 


The term diag(pj indicates a diagonal matrix, with all diagonal terms equal to the 
density of the k-th ply. The kinetic energy in eq. (12) includes contributions of rotational 
inertia. As seen, the proposed laminate damping mechanics are general and they can 
handle any laminate configuration in terms of composite plies, pfy angles, laminate 
stacking sequence, and thickness. 
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33. Simply-Supported Composite Plates 


Based on the previously described laminate mechanics, exact solutions of the damped 
dynamic characteristics for specialty composite plates have been obtained and 
summarized herein for the sake of completeness^^ These exact and computationally 
inexpensive predictions of static and d)oiamic characteristics (modal damping, natural 
frequencies) provided valuable insight in the mechanics of the problem and establish a 
baseline for validating approximate numerical methods. Most other laminations, boundary 
conditions, and structural configurations require approximate discretized solutions, and 
such finite element based computational mechanics are described in the next section. 

For a rectangular a by /3 simply supported (SS) composite plate with negligible coupling 
= A26=0, = B^26=0, = E>*“2<;=0), the following Navier fundamental 

solutions form a complete set of mode shapes in the x-y plane: 

= tCcos(ox)sin(fey)e'"' 

= C,sin(^)cos(6y)e'“' 

= PCsin(ax)sin(i>y)e'- 
= Uj^cos(ax)sm(byy'^‘ 

Where, a= mn/a, b= mr/jS, and the fundamental solutions in eq. (14) are chosen because 
they satisfy the boundary conditions of the simply-supported plate. 

Combination of eqs. (4,5,14) yields the modal mid-plane and generalized strains as 
separable functions of x, y coordinates, time, and amplitudes. 


where, the terms in matrices [B^ J are sinusoidal functions of x, y coordinates and mode 
order. The amplitude displacement vectors are {IT} = {U° , V® , W®}^ and {U*} = {U\ 

v>r. 


By substituting eqs. (15) into the strain and kinetic energy expressions for the laminate, 
eqs. (10) and (12) respectively, integrating over the plate area, and applying Lagrangian 
dynamics, the modal analysis (free vibration) solution of the plate takes the form: 
where, the through-the-thickness modal displacements are {U^> = {U®^; 
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( 16 ) 


and subscripts m, n indicate the mode order. Numerical solution of this eigenvalue 
problem provides the natural frequencies o)nj„ and the through-the-thickness modes {U^} 
for each order mn of plane modes, in the context of eqs. (14). 

The modal damping associated with the mn-th vibration mode rjnm 




2w 



where AS, „„ and S | are respectively the dissipated and maximum laminate strain 
energies, respectively of the associated mode. 
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3.4. Finite Element Method 


Computational finite-element based mechanics for the damping analj^is of thick 
composite structures of general laminations and shapes were developed, based on the 
DLDT and are summarized herein. In view of the generalized laminate strains and 
stresses described in the previous sections, the variational formulation of the equilibrium 
equations for the candidate structure takes the form: 

N ^ 

i-i ( 18 ) 

f; £ ")<M .Siii’f,./ 

/-I W*1 

For viscoelastic constituents, the constitutive viscoelastic law between generalized stresses 
and strains at the laminate level may be eventually expressed by the following convolute 
products. 


N 

N\t) = [i(t)]*e^(f)+S [B'"(f)]*e"(0 

N 

«*1 

where * indicates a convolution operation, for example. 


(19) 


Hence, the first left-hand integral in eq. (18) represents both stored and dissipated strain 
energies during the virtual displacement. The elastic laminate relations will result from 
eq. (19) if time effects are neglected. Also, viscous damping is a special case of eqs. (18). 
In the frequency domain (harmonic vibration), eqs. (19) will result in equivalent 
expressions of complex laminate stiffness, with storage and loss stiffiiess matrices 
described in eqs. (8) and (11) respectively. 


The variational relationship may be discretized to an approximate equivalent system of 
dynamic equations, if the reference and generalized displacements are discretized as 
follows: 


where, subscript i indicates the i-th node, and 4 >, are interpolating functions. The 
reference and generalized strain shape functions, [R**]; and [R]j interpolating the strains 
to the nodal displacements, are obtained by differentiating eq. (21) in accordance with 
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( 21 ) 


nodes 

i-1 
nodes 

i*l 

eqs. (4,5). Then the damping or loss element matrix [C^jj] for nodes i and j takes the 
form: 




Khsi]Rj ... R^W^Rj 

R^^[D!]RJ ... R,W}Rj 

• • * • 

• « « 

••• R^i^r^Ri 




( 22 ) 


where is the area of the element. The area integration is performed numerically The 
stiffness and mass matrices of the element, [K^jj] and [M^j] respectively, have analogous 
forms. Subsequently, the dissipated and maximum stored strain energies of the element 
are given by: 


AS,, =1/2 


S„.=l/2 


(23) 


and the modal damping of the n-th mode is expressed as the ratio of the dissipated 
over the maximum modal strain energy of the structure. 


1 E AS . 
■ 27C E S. 

. eijn 


(24) 


In this context, a bilinear four-node plate element was developed. Selective integration 
was utilized to avoid overstiffening of the element in the case of low thickness. Such 
overstiffening was observed with full integration. 
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LAMINATE 



THROUGH-THE-'miCKNESS 

displacements 



CLASSICAL LAMINATE 
PLATE THEORY 
ffHIN LAMINATES) 



DISCRETE LAMINATE 
PLATE THEORY 
(THICK LAMINATES) 


Fig. 3.1 (a) Unidirectional composite, (b) Kinematic assumptions of laminate 

damping theories; classical laminate damping theory (CLDT) and discrete 
layer damping theory (DLDT). 
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4. DAMPING PREDICTIONS IN THICK COMPOSITE STRUCTURES 


To demonstrate the quality and versatility of the composite mechanics, applications on 
Gr/Epoxy composite plates and beams were performed. When applicable, numerical 
results obtained with the aforementioned element were compared with the exact 
solutions. In all other cases the predictions were compared with results obtained using a 
triangular plate element based on classical laminate damping theory (CLDT)“. Both 
elements will be referred to respectively as DLDT and CLDT finite elements thereafter. 
The composite material used in all applications was HM-S Graphite/Epoxy of 50% fiber 
volume ratio (FVR). Measured data for this composite are provided in Ref. 9. The 
properties of the the fibers and the matrix were backcalculated using micromechanics’ to 
match the measured properties of the composite. Unless otherwise stated, the nominal 
thickness of each composite ply was taken as 0.254 mm (0.01 inches). 

4.1. [ 07904 ], Simply-Supported Composite Plate 

Fig. 4.1 shows predicted natural frequencies and modal damping values of the 
fundamental mode of a [O 79 O 4 ], square composite SS plate as functions of thickness 
aspect ratio (a/h), where t is the plate thickness. The modal characteristics of the first 4 
modes, that is, modes (1,1), (1,2), (2,1), and (2,2) are shown respectively. It is pointed 
out, that this laminate has different flexural damping and stiffness characteristics in x and 
y directions, thus, modes (1,2) and (2,1) have different natural frequencies and modal loss 
factors. The finite element model consisted of a uniform mesh of 8 by 8 element 
subdivisions. Qearly, there is excellent agreement between the exact DLDT solutions and 
the finite element model in the shown thickness regime. The exact CLDT solution is also 
shown, to illustrate the significant differences between the two theories in damping 
prediction at high thicknesses. The differences between DLDT and CLDT predictions 
provide the estimate of interlaminar shear damping. Finally, the results calculated using 
the DLDT finite element converge to the exact CLDT predictions as the plate thickness 
is decreased. 

Fig. 4.2 shows the effect of temperature on the fundamental modal damping and natural 
frequency of [O 79 O 4 ], and [0/90]4,, 203 mm by 203 mm (8 in by 8 in) square plates. 
Predictions of the fundamental SDC and natural frequency of each plate, subjected to 
uniform through-the-thickness temperature variations are plotted as functions of 
temperature. Interestingly, the CLDT significantly underestimates the fundamental SDC 
at higher temperatures. The underestimation is slightly lower in the second plate. It is 
recalled, that in higher temperatures both the shear modulus of the 90-deg. plies 
decreases and the shear damping increases respectively. The temperature variation has 
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little effect on the natural frequency. 


4.2 [0/90/45/-45]5 Free-Free Composite Beam 

Fig. 4.3 presents the modal damping and natural frequencies of the first and second 
bending modes of an unsupported (free-free) [0/90+ 0/45+ 0/-45+0]s thin beam 
(l/h=125). This laminate configuration was used previously in both experimental and 
analytical studies®’®. The CLDT model has failed to capture the portion of damping due 
to interlaminar shear stresses and relative rotation between adjacent plies. Hence, this 
case serves as a demonstration of the dependence of damping on the specific laminate 
lay-up only, and also as an example of the capacity of DLDT to capture both flexural 
and shear damping. 

Equivalent element meshes consisting of 15 by 4 element subdivisions were used to 
predict the dynamic characteristics of the beam. The dimensions of the beam were 
approximately equal to specimens used by Ni and Adams®, that is, 200 mm (8 in) long, 12 
mm (0.5 in) wide, and 1.60 mm (0.064 in) thick. As seen in Fig. 4.3, the DLDT damping 
predictions are higher than the CLDT prediction and correlates very well with the 
reported measurements®, hence, the new theory has succeeded in capturing the 
interlaminar shear damping contributions in the range of high interlaminar shear stresses 
and relative rotation between outer plies (0- 45°, 135°). It is emphasized, that in this case 
the interlaminar shear damping is mostly the result of inhomogeneities in the material 
anisotropy between the plies of alternating fiber orientation; therefore, the case study 
serves as another example of the superiority of the method. 

43. [ 04 /- 04 ls Simply-Supported Composite Plate 

The lower modal damping values and natural frequencies predicted for a simply 
supported [ 04 /- 04 ]s GryEpoxy plate are shown in Fig. 4.4. A uniform 8 by 8 finite element 
mesh was again used to model the plate. The modes are identified with their mode shape 
at 0=0. The plate is rather thin (a/h=100), therefore, the predictions with both DLDT 
and CLDT finite elements were expected to be similar. Yet, the DLDT element resulted 
in higher predictions of damping with respect to the CLDT element, as a result of the 
additional shear contributions. The difference in damping predictions, which indicates the 
contributions of shear damping, increased with the order of the mode. It seems that the 
interlaminar shear damping contributions are slightly higher for fiber orientations around 
±45°, which corresponds to lay-ups with higher interlaminar shear stresses, indicating the 
advantage of DLDT in representing interlaminar shear damping. 
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To investigate further the dependence of interlaminar shear damping on the laminate 
lay-up, the modal damping of the previous 3 modes is shown in Fig. 4.5 for plates of 
various thickness. The damping predictions of the CLDT model, which are independent 
of thickness, are also plotted. The corresponding natural frequencies are shown in Fig. 
4.6. Keeping in mind that the difference between DLDT and CLDT results indicate the 
magnitude of shear damping, there is defimte coupling between lammate lay-up and 
thickness, especially in the second and third mode. The interlaminar shear damping is 
again higher at ±45 degs. where the interlaminar stresses are highest. In this manner, the 
results provide additional confidence and value to the damping mechanics. 


4.4. [BfiZ-Osls Cantilever Composite Beams 

Figs. 4.7 to 4.9 present the predicted fundamental modal damping and natural frequency 
of [BsZ-Bsls cantilever beams of thickness aspect ratios (1/h) 60, 30, and 15 respectively. 
The width and thickness of the beams were 25.4 mm (1 in) and 0.508 mm (0.2 in) 
respectively. The same uniform mesh consisting of 15 by 4 element subdivisions was used 
in all cases. For ply angles between 0 and 60 degs., the DLDT model 3 delded significantly 
higher predictions of damping and lower predictions of natural frequency in the cases of 
thicker beams (1/h =15). These differences were attributed to the high interlaminar 
stresses at the higher thickness beams. The coupling between lay-up and thickness of the 
beam is also illustrated in Fig. 4.9. The results demonstrate again the quality of the 
method in predicting shear effects due to thickness, laminate lay-up, or combinations of 
both. 
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Fig. 4.1 Modal characteristics of a [04/904]^ Gr/Epojqr simply-supported plate. 
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Fig. 43 Modal characteristics of a [0/90+ 0/45 +0/-45+e]s cantilever beam (l/h=125; 
1=8 in). First and second bending modes. 
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Fig. 4.9 Fundamental modal characteristics of a Gr/Epoxy cantilever beam. 

(l/h= 15) 
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5. COMPOSITE PLATES WITH EMBEDDED COMPLIANT DAMPING LAYERS 


5.1. Experimental Correlations 

To validate the accuracy of the method and experimentally verify the effects of damping 
layers on the dynamic response of composite plates, a composite plate with interlaminar 
damping was fabricated in-house and tested in-house^ (see Fig. 5.1). The plate was 11 in 
by 11 in (28 cm by 28 cm) and consisted of plies of 0.6 FVR T300/934 Graphite/Epo)q^, 
combined with layers of ScotchDamp ISDllO damping film (3M Corp., St. Paul, MN), laid 
up in a [02/902/i/902/02]s configuration, where i denotes the damping layer. The shear loss 
factor and moduli of the damping polymer, as approximated with eqs. (24) and (25), are 
shown in Fig. 5.2 for the frequency range of interest, and represent average values of the 
polymer at 30® C provided in the manufacturer’s data sheet. The nominal thicknesses of the 
Gr^p tape and damping layers were 0.006 in (0.15 mm) and 0.005 in (0.13 mm), 
respectively. The plate was laid up by hand, co-cured at 350° F (175° C) and 50 psi (345 
KPa), and finally cut to size. 

Frequency response functions for the composite plate were obtained fi-om impulse tests. A 
special fijrture with “knife-edge” supports was designed to support as closely as possible the 
plate in the SS configuration described by eqs. (16), and to minimize the friction forces. The 
plate was impacted at the quarter point (x/a, y//3)= (0.25, 0.25) with an instrumented 
hammer, and the acceleration was measured with a miniature accelerometer (1 gm mass) 
at points (x/a, y/^)= (0.5, 0.5) and (x/a, y//3)= (0.25, 0.75). The signals from the hammer and 
the accelerometer were recorded with a high speed digital data acquisition system, and was 
then processed using FFT software to obtain the frequency response functions of the plate. 
The accelerance of the plate at the previously mentioned points is shown in Fig. 5.3. 

The predicted response was synthesized from the calculated modal damping and natural 
frequencies using a modal superposition technique (Fig. 5.3) with frequency effects included 
as described in the previous section. With the exception of the first mode, the predicted 
response reproduces sufficiently the measured one. The mismatch and higher damping of 
the first mode was also observed in the test of a monolithic aluminum plate, hence, it was 
attributed to imperfections in the experimental supports. There is some underestimation of 
the damping in the predicted response, which may be attributed to firiction and air damping 
present in the experiment. However, both predicted and measured results fall within the 


^ The author wants to thank Mr. Peter Addente for fabricating the composite plates 
and Dr. Michael J. Pereira for testing the plates. Their contributions are 
gratefully acknowledged. 
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range of uncertainty predicted front the upper and lower bounds of scatter in the 
manufacturer’s data for this polymer film damping. The agreement is improved at higher 
modes, because the effects of the supports are less dominant compared to the shear 
damping contributions of the damping material. 


52. Parametric Studies 

This section presents applications of the developed methodology on square laminated 
Gr/epoxy simply-supported (SS) plates with interlaminar damping layers. The dimensions of 
the plate were assumed 11.8 in by 11.8 in (30 cm by 30 cm). The composite material was 
HM-S graphite/epojq^ with elastic and damping properties given in ref. 6. Unless otherwise 
stated, the fiber volume ratio (FVR) of the composite plies were 0.50. The shear modulus 
and shear loss factor of the interlaminar damping layer were 0.9 kpsi (6.2 MPa) and 1.02 
respectively, and represent average values of a typical commercial damping polymer 
(ScotchDamp ISD112, 3M Corp)^^ at room temperature. Experimental data regarding the 
normal loss factor and Young’s Modulus of the damping polymer was not available. Because 
the present theory includes provisions for normal damping and modulus, a Poisson s ratio 
equal to 0.3 and a very low normal loss factor (0.016) were assumed for the damping 
polymer to ensure minimal contribution on the damping of the plate. This assumption was 
further reinforced by parametric studies which quantified the insensitivity of the modal 
characteristics (natural frequencies and damping) of the plate to the normal damping and 
the Poisson’s ratio of the damping layer. The effects of frequency on the properties of the 
damping polymer were neglected here, in order to isolate the effects of composite anisotropy 
and inhomogeneity on the performance of the damping layers. The nominal thickness tj of 
the damping layer was equal to the typical thickness of a composite ply t, = 0.005 in (0.13 
mm). In all cases a concentrated vertical load of 100 lbs (445 N) was applied at the center 
of the plate. Linear interpolation functions F’(z) were used to calculate the in-plane 
displacements through-the-thickness of the laminate. 

The ply-stacking sequence is represented using the traditional convention, but interlaminar 
damping layers of thickness tj are identified with symbol i. For comparison purposes, 
calculated results for aluminum plates are also presented, in such case, an aluminum layer 
of thickness equal to the thickness of k plies will be identified as Al^. All laminate 
configurations were symmetric. The modal loss factor is normalized by the shear loss factor 
of the polymer, and will be referred in the next paragraphs as normalized modal damping. 

Predicted modal damping, natural frequencies, and static deflections at the center of an 
aluminum (Al 4 /i/Al 4 )s and various Gr/epoxy damped plates, ( 04 /i/ 04 )s, (04/i/904)s and 
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(02/902/1/02/902)5 are shown in Table 5 . 1 . All plates have high modal damping values. The 
laminates with the higher anisotropy in the composite sublaminates, ie. (04/i/04)s and 
(04^904)8, exhibit the higher damping in the fundamental mode. The (04/1/904^ laminate 
exhibits the higher anisotropy variation through-the-thickness among all plates, and this 
seems to have resulted in some higher modal damping values. Interestingly, this laminate has 
comparable static stif&iess and natural frequency with the ( 02 / 902 /i/ 02 / 902 )s laminate, which 
illustrates the possibility for higher damping values in some modes without significantly 
penalizing other structural characteristics. The high damping of the unidirectional composite 
laminate was attributed in part to the lower stiffness of the plate. The composite plates have 
higher static deflections but lower areal density than the aluminum plate, as a result, they 
exhibit slightly higher fundamental frequencies. The damping of the higher modes varies 
widely and depends on the respective mode shape. The variation is more predominant for 
the composite plates. 

The variations of the first modal damping, natural firequenqr, and static deflection at the 
center of the plate are plotted in Fig. 5.4 as functions of the damping layer thickness. The 
thickness of the composite sublaminates was kept constant, therefore, the total laminate 
thickness was varied. The layer thickness tj is normalized by the nominal ply thickness tj= 
0.005 in. For layer thicknesses in the range of practical interest, the laminate with the higher 
anisotropy through-the-thickness (04/i/904)s provides consistently higher damping than the 
aluminum (Al4/i/Al4)s and ( 02 / 902 /i/ 02 / 902 )s laminates. This results are consistent with the 
expectation that composite laminates will induce higher interlaminar shear stresses in the 
interlaminar layer, in order to balance the elastic anisotropy and nonuniformity of adjacent 
plies, as opposed to traditional isotropic laminates where interlaminar stresses will only 
balance the bending stresses through-the-thickness. For all laminates, the modal damping 
increased with thickness and tends to reach a maximum at higher thickness values. It was 
observed that for plates of lower aspect ratio alh the peak in the modal damping shifts 
towards lower thicknesses. Finally, it seems that a critical damping layer thickness exists (in 
this case approximately equal to the ply thickness), such that, damping layers with less or 
equal thickness will add low but significant damping with negligible effects on the structural 
properties (static deflection, natural frequency). 

Fig. 5.5 shows the predicted variations of modal damping, natural frequency and static 
deflection with the aspect-ratio of the plate (a/h). As was expected, the modal damping of 
all plates increases at lower aspect ratios because shear prevails over flexure. Interestingly, 
the damping reaches a peak at very low aspect ratios and then decreases again. Intuitively, 
at very low aspect ratios through-the-thickness deformations will gradually prevail over shear 
and the effectiveness of the damping layer will be decreased. Nevertheless, the kinematic 
assumptions of the DLDT are not expected to be valid at the regime of very low aspect 
ratios because the present approach assumes negligible through-the-thickness strains. Both 
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natural frequency and static stiffness consistently increase with lower aspect ratios. 

Fig. 5.6 presents the variation of modal damping, natural frequency, and static deflection as 
functions of the FVR of the inner and outer composite sublaminates with respect to the 
damping layers. Fig. 5.6 demonstrates one unique advantage of composite laminates, that 
is, the additional tailoring capacity they can introduce. In both laminates, the FVR variation 
of the outer sublaminate has a definite effect on the modal damping of the plate. The 
damping increases almost linearly with the FVR, indicating that the most important 
composite parameter is the longitudinal modulus. Since both natural frequency and static 
stiffness also increase with the FVR of the outer sublaminate, the obvious choice is to have 
outer sublaminates of very high FVR. The modal damping is less sensitive to FVR variations 
of the iimer sublaminate, and decreases slightly with FVR. 
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Table 5.1. Static and dynamic characteristics of composite plates with interlaminar layers. 
(tj=0.005 in; numbers in parenthesis indicate mode order) 



\ 1 

(Al4/i/Al4)s 

(OyWs 

(04/i/904)s 

(02/902/i/02/902)s 

Areal Density, 10'^ Ib/in^ 





7.95 

4.92 

4.92 

4.92 

Static Deflection (Center), in 




0.334 

0.604 

0.496 

0.456 

Normalized Modal Damping 



Mode 1 

0.131 (1,1) 

0.153 (1,1) 

0.153 (1,1) 

0.098 (1,1) 

2 

0.259 (1,2) 

0.126 (1,2) 

0.066 (1,2) 

0.266 (1,2) 

3 

0.259 (2,1) 

0.104 (1,3) 

0.433 (2,1) 

0.252 (2,1) 

4 

0.346 (2,2) 

0.110 (1,4) 

0.355 (2,2) 

0.277 (2,2) 

5 

0.374 (1,3) 

0.370 (2,1) 

0.033 (1,3) 

0.438 (1,3) 

6 

0.374 (3,1) 

0.355 (2,2) 

0.218 (2,3) 

0.395 (2,3) 

7 

0.415 (2,3) 

0.319 (2,3) 

0.598 (3,1) 

0.381 (3,1) 

8 

0.415 (3,2) 

0.134 (1,5) 

0.553 (3,2) 

0.380 (2,3) 

9 

0.448 (1,4) 

0.277 (2,4) 

0.028 (1,4) 

0.545 (1,4) 

10 

0.448 (4,1) 

0.249 (2,5) 

0.438 (3,3) 

0.418 (3,3) 

Natural Frequencies, Hz 




Mode 1 

115.7 

119.9 

119.9 

124.0 

2 

265.0 

152.7 

211.8 

272.6 

3 

265.0 

230.1 

354.8 

357.7 

4 

395.2 

351.2 

402.4 

437.7 

5 

474.5 

383.4 

409.3 

495.2 

6 

474.5 

402.3 

541.7 

606.7 

7 

585.6 

447.7 

635.2 

688.7 

8 

585.6 

508.7 

667.9 

736.7 

9 

722.5 

531.4 

697.6 

747.2 

10 

722.5 

656.5 

766.7 

851.6 
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Fig. 5.1 Experimental configuration showing composite plate in knife-edge supports. 
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Fig. 5.2 




Mechanical properties of the ISDllO damping polymer at 30° C. Shear 
modulus and shear loss factor. 
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Fig. 53 Frequency response of a (02/902/i/902/0^s T300/934 plate with ISDllO 
damping layers. 
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Fig. 5.4 Effect of damping layer thickness. Alum.= (Al 4 /i/Al 4 )s; Lam. 1- (04/i/904)s; 
Lam. 2= 


36 








6 . DAMPING ENHANCEMENT IN SPECIALTY COMPOSITE STRUCTURES 


The developed structural damping mechanics were also applied on the analysis of 
damping in various composite plates with embedded compliant interlaminar damping 
layers. The objectives of the case studies were: (1) to further evaluate the quality of 
method by correlating predictions with both analytical and experimental results; and ( 2 ) 
to explore and quantify unique advantages in improving the vibroacoustic characteristics 
of composite structures by embedding compliant layers between plies with strong 
interlaminar inhomogeneity. These applications demonstrate the potential to 
engineer/tailor such specialty composite laminates for high dynamic performance. 

The selected materials were 0.60 FVR T300/934 composite, and ScotchDamp ISDllO 
damping film (3M Corp., St. Paul, MN) as interlaminar layers. The shear modulus and 
shear loss factor of the damping polymer were assumed 5.5 kpsi and 0.16 respectively, 
and they represent average values of the polymer at 500 Hz and 30® C as provided by 
the manufacturer’s data sheet. Due to non-existing data, the Young’s modulus of the 
damping polymer was assumed 15.4 kpsi and the normal loss factor was assumed to be 
negligible (0.016). Parametric studies have indicated strong insensitivity of modal damping 
and frequency to the dilatational properties of the damping polymer. The damping 
properties of the 934 epoxy were assumed equal to typical damping values of epoxy 
matrices. Unless otherwise stated, the nominal thicknesses of the Gr/Ep plies and 
damping layers were 0.0058 in (0.15 mm) and 0.005 in (0.13 mm), respectively. 


6.1. [ 04 / 1/90418 Simply Supported Plate. 

The quality and accuracy of the developed structural damping mechanics were first 
investigated by comparisons with exact analytical solutions. Fig. 6.1 show^ the modal 
damping and natural frequencies of a [ 04 / 1 / 904 ]$ simply supported plate predicted by the 
discrete finite element method and the exact plate solutions developed based on the 
discrete laminate damping theory [10]. The finite element analysis involved an 8 by 8 
uniform mesh subdivision. As seen in Fig. 6.1, there is excellent agreement between the 
finite element based and exact damping predictions. The accuracy of the predictions is 
slightly but gradually reduced at higher modes, because of the limitation of the 
implemented mesh in representing the higher mode shapes. A more refined mesh would 
readily improve the predictions at higher modes. 


6.2. Angle-ply Unsupported Plates. 
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The predicted dynamic characteristics of two symmetric angle-ply plates with compliant 
damping layers are also shown in Figs. 6.2 and 6.3. The laminate configurations were 
[452/-452/i/452/-452], and [22.52/-22.52/i/22.52/-22.52]3 respectively. In-house experimentaUy 
measured damping values and natural frequencies for these plates at a free-hanging 
configuration are also shown in these figures. The composite plates were fabricated by 
hand lay-up, co-cured at 350° F (175° C) and 50 psi (345 KPa), and finally cut to 28 cm 
by 28 cm (11 in by 11 in). 

In both cases, the correlation between measured and predicted modal damping and 
corresponding natural frequencies appears very good, especially in the low frequency 
regime. Predicted and measured results for composite plates of identical laminations but 
w/o damping layers are also shown in Figs. 6.2 and 6.3. It is pointed out, that the sole 
purpose of plotting these results is to compare and demonstrate the obtained 
improvements in the damping characteristics of the various plates, rather than comparing 
the accuracy of the method, since damping data for the composite system were not 
available at the time of this paper. 


63. Angle-Ply Laminates 

The predicted loss factors and corresponding natural frequencies of the first five modes 
of two simply-supported plates, with laminations [02/-02^®2^‘®2l* [02/-02/^‘®2/‘^®2]s 

respectively, are shown in Figs. 6.4-6.8 as functions of the ply an^e 0. The modal loss 
factors and natural frequencies of plates with identical geometrical and laminate 
configurations, but without damping layers (undamped), that is [02/-02/02/-®2]s l®2/* 

02/-02/+02]s respectively, are also plotted in the same Figures. 

Qearly, the addition of the damping layers drastically improved the damping capacity of 
the composite plates. But most importantly, the predicted results demonstrate that 
increased through-the-thickness inhomogeneity and anisotropy in the composite 
sublaminates may significantly improve the damping of the composite structure. This 
damping enhancement is further signified by the observation that the predicted damping 
of the composite plates without damping layers is typically reduced for ply angles 
between 0 and 90°. Finally, the damped plates exhibit similar or higher natural 
frequencies than the undamped plates between 0 and 90 degs, hence, this enhancement 
of damping does not seem to reduce the natural frequencies and stiffiiess. 

As seen in Figs. 6.4-6.8, the improvements in the damping are increased at higher modes, 
for example mode 1 (1,1), mode 2 (1,2), mode 3 (1,3), and mode 5 (2,2). The numbers in 
parenthesis indicate the mode order at 0=0 degs. The damping increases, indicate the 
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interrelation between laminate configuration and mode shape. A significant part of the 
damping improvements were induced by altering the anisotropy in the restraining 
composite sublaminates. It seems though, that alternating fiber orientations in the plies 
which contain the damping layer produce some additional damping improvement. It is 
also interesting to note, that as the stiffness of the restraining sublaminate remains the 
most critical mechanism for enhancing shear damping in the interlaminar layer, not all 
modes will exhibit the previous damping improvements. This is illustrated in mode 4, 
which corresponds to mode shape (2,1) and exhibits its highest damping at 6=0 degs, 
that is, when the outer restraining sublaminates have the higher flexural stiffiiess in the 
direction of the flexure. 

The presented results, illustrated the accuracy and versatility of the structural damping 
mechanics in predicting the damped dynamic characteristics of specialty composite 
structures. The applications have also demonstrated the potential to "engineer" composite 
laminates with embedded interlaminar damping layers to provide higher and tailorable 
damping than traditional isotropic materials. The results finally demonstrate the 
opportunity to tailor and "engineer" composite structures for improved dynamic 
characteristics in terms of damping and natural frequencies. 
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Fig. 6.1 Modal characteristics of a (04/i/904)s T300/934 simply supported plate with 
ISDllO damping layers. 
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Fig. 6.2 Modal characteristics of a [452/-452/i/452/-452]s T300/934 simply supported 
plate with ISDllO damping layers. 
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Fig. 63 Modal characteristics of a [22.52/-22.52/i^2.52/-22.52]s T300/934 simply 
supported plate with ISDllO damping layers. 
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Fig. 6.4 Modal loss factor and natural frequency of angle-ply T300/934 simply 
supported plates; mode 1. 
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Fig. 6.5 Modal loss factor and natural frequency of angle-ply T300/934 simply 
supported plates; mode 2. 
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Fig. 6.7 Modal loss factor and natural frequency of angle-ply T300/934 simply 
supported plates; mode 4, 
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Fig. 6.8 Modal loss factor and natural frequenqr of angle-ply T300/934 simply 
supported plates; mode 5. 
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7. OPTIMAL DESIGN OF THICK COMPOSITE PLATES 


7.1 Assumptions 

This section presents preliminary evaluation of the developed optimal design 
methodology for composite plates. For this purpose, the optimal design of a thick 
composite plate was investigated as a minimal weight problem. Although the design of 
composite structures for optimal dynamic performance is really a multi-objective task, the 
present initial investigation maintains simplicity in the objective function. Increases in 
composite damping may typically result in stiffness/strength reductions and/or mass 
addition, for this reason, the minimization of weight subject to deflection, natural 
frequency and modal damping constrains is a reasonable first formulation. 

A constrained single-objective problem is described in the following mathematical form: 

min {Fj(d)} (25) 

subject to lower and upper bounds on the design vector d and inequality constraints 
G(d): 

{d^} ^ {d\ ^ [d^] 


{G(d)] ^ 0 

In the rest of the paper, upper and lower values are represented by superscripts L emd U 
respectively. 

The objective function in eq. (9) is assumed to be the specific weight per unit area of the 
plate. The design vector may include fiber volume ratios (FVRs), ply angles, and 
thicknesses of various sublaminates. Performance constraints include upper bounds on 
static deflections u®, 

{«*} 5 

lower and/or upper bounds on natural frequencies {!„,„}> 


(26) 

(27) 
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if-) i ^ if”) 


(29) 


where, m, n indicate mode order. Additional lower bounds are imposed on each modal 
damping of the composite plate to ensure that the passive damping of specific modes 
remains or is "pushed" above desirable limits. 

{ii„) 1 (<i‘i 

Constrains on static and dynamic stresses of composite plies are not included in the 
present study, but their inclusion is planned as future work. 

The constrained optimization mentioned above are solved with the modified feasible 
directions non-linear programming method^*. 


12 Optimal Design of a [0/90] ^ Plate 


Evaluations of the method on Gr/Epoxy cross-ply [0/90]s composite plates were per- 
formed. The fiber volume ratio (FVR) and the thickness of the 0 and 90 degs 
sublaminates were optimized. The initial FVRs were assumed equal to 0.50. The initial 
thickness of each sublaminate was assumed equal to 0.040 in, that is, the initial thickness 
of the plate was h=0.16 in. The plates where assumed square (R=a/j8=l), but various 
length values (a) were considered, in order to investigate the effect of thickness-aspect 
ratio (ct/h) on the resultant optimal designs. An out-of-plane concentrated load of 100 lbs 
at the center of the plate was applied in all cases. 


Apparently, as various values of the plate span a are considered, the natural fi’equencies 
and the static deflections will also change, hence, non-dimensional expressions for these 
expressions were utilized. In the context of CLPT solution^^, for constant R, the non-di- 
mensional expressions of transverse displacement, and natural frequencies may be 
defined as: 


H<x,y) = 


w(Jg>y) 

Aa^ln* 





2|«2 


n la 


(31) 


These expressions were utilized to scale the bounds on constraints (12) and (13) when 
different values of a were considered, such that geometric similarity in the context of 
CLPT was maintained. The objective function and modal damping were already non- 
dimensional. As an example, for a plate of a=/3=10 in, an upper bound on the static 
defection at the center of the plate was set equal to 0.050 in; lower bounds of 300, 500, 
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1000, 1000, and 1000 Hz were imposed on the first five natural frequencies of modes 
(1,1), (1,2), (1,3) and (2,2); and lower bounds of 0.24% was imposed on the 
corresponding loss factors of the previous five modes. For other values of span length, 
the constrains were scaled as previously described. 

Figs. 7. 1-7.5 illustrate the effect of thickness aspect ratio on the resultant objective 
function and corresponding optimal design parameters. The corresponding a/h was 
calculated at the initial designs. Because of the scaled objective and constrains, the 
optimal designs based on CLPT appear independent to the thickness aspect ratio. 

Clearly, the resultant optimal designs based on DLPT are far more sensitive to aspect 
ratio than the actual vibrational characteristics (Fig. 4.1). Only at very high aspect ratios 
(a/h « 200) comparable optimal designs were attained. 

The results illustrate, that at least for the current application, the optimal design of 
composite plates is very sensitive to the composite mechanics used. Design methods 
based on classical composite mechanics may produce significant errors even in case of 
intermediate to low thicknesses. Moreover, the results demonstrate the quality and 
superiority of the present design method which is based on a more accurate laminate the- 
ory. Additional evaluations are planed to be performed in the near future. 
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Fig. 7.4 Optimal thickness of 0° sublaminate for various initial thickness aspect 
ratios (a/h). 



Fig. 7.5 Optimal thickness of 90° sublaminate for various initial thickness aspect 
ratios (a/h). 
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CONCLUSIONS AND FUTURE WORK 


The development of novel unified damping mechanics for the accurate analysis and 
design of composites structures was summarized. The work emphasized the need of 
damping mechanics applicable to represent not only thin, but also thick composite 
structures, and specialty composite structures with embedded compliant interlaminar 
layers. The effects of elevated temperature in the (tynamic response of the composite 
structure were also included. The more significant results, accomplishments and 
prototype software developed as part of this effort were enumerated. 

Numerous applications and evaluations of the mechanics illustrated their merit and 
unique advantages in the case of higher thicknesses, elevated temperatures and 
interlaminar inhomogeneities. The importance of the mechanics in the optimal design of 
composite structures was also demonstrated. 

Future work is planned in the following areas: 

(1) development of a three-dimensional damping theory is planned by 
incorporating a variable through-the-thickness displacement fields. This theoiy will 
establish the upper bound of accuracy, and also will entail the capacity to model 
the interlaminar normal stress, local impact phenomena, and through-the-thickness 
wave propagation (vis. noise transmission and acoustic characterization); 

(2) the inherent capacity of the theoiy to represent interlaminar damage in the 
material will be explored and the effects of delamination on the damping and the 
dynamic characteristics will be quantified. 

(3) extension of the mechanics to model health monitoring conditions in "smart" 
composite structures with embedded sensors in the form of layers and/or fibers. 
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APPENDIX I: Nomenclature 


[A], m pf“] 
[A,], [BA PA] 

[A„], PA PJ“] 
[C] 

E 

P(z) 

G 

h 

[K] 

1 

[M] 

N°, Nj 

[Q] 

PI 

s 

U, V, w 


x,y, z 
a, jS 
AS 
{€} 
n 


p 

0 


ilr 


Extensional, coupling and flexural/shear laminate stiffiiess matrices. 

Extensional, coupling and flexural/shear laminate damping matrices. 

Extensional, coupling and flexural/shear laminate mass matrices. 

Damping matrix, structural 

Normal modulus 

Interpolation function 

Shear modulus 

cross-section thickness 

Stif&iess matrix, structural 

beam length 

Mass matrix, structural 

Midplane and generalized stress resultants 

Composite stiffness matrix 

Strain shape function 

Maximum stored strain ener^, per cycle 

Axial, transverse and through-the-thickness displacements 

respectively 

Axial, transverse and through-the-thickness structural axes. 

Plate dimensions along x and y axes, respectively 

Dissipated strain energy, per cycle 

Engineering strain tensor 

Loss factor 

Poisson’s ratio 

Mass density 

Shape function 

Natural frequency 

Specific damping capacity 


Subscripts 

c Composite (structural coordinates) 

e Element 

1 Composite (material coordinates) 

L Laminate 

m Mode order, x direction 

n Mode order, y direction 

1 Normal longitudinal (direction 11) 

2 Normal transverse, in plane (direction 22) 

3 Normal transverse, out of plane (direction 33) 
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4 

5 

6 


Interlaminar shear (direction 23) 
Interlaminar shear (direction 13) 
In-plane shear (direction 12) 


Superscripts 

o 

j 


Midplane 

Generalized 
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APPENDIX II: Fundamental Composite Mechanics 


Off-axis ply stiffness matrix [QJ: 


[<?j= 


where [QJ is the stiffness matrix in the material axes. The transformation matrices are: 


[R] 


m2 

n2 

0 

0 

0 

2mn 


m2 

0 

0 

0 

-2mn 

0 

0 

1 

0 

0 

0 

0 

0 

0 

m 

-n 

0 

0 

0 

0 

n 

m 

0 

-mn 

mn 

0 

0 

0 

m2-n' 


m = cos6, n = sin0 


[R]-^ = [i?(-6)] 

On-axis damping matrix: 

[hj “ [ ^llV ^12Z> ^133^ ^144'> ^155> ^166 1 

The off-axis ply damping matrix is: 

[iij= [mtijEiq'' 
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^cJI 

^cl2 

0 

0 

0 

^cl6 

%2I 

^c22 

0 

0 

0 

%26 

0 

0 

^c33 

0 

0 

0 

0 

0 

0 

Vc44 


0 

0 

0 

0 

%54 

%S5 

0 

%61 

^c62 

0 

0 

0 

'^c66. 
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